Constraints on the Scalar- Tensor theories of gravitation 

Primordial Nucleosynthesis 



A. Serna and J. M. Alimi 

DAEC, Ohservatoire de Paris-Meudon, 92195-Meudon, France 
(October 31, 1995) 

Abstract 

We present a detailed calculation of the light element production in the 
framework of Scalar- Tensor theories of gravitation. The coupling function uj 
has been described by an appropriate form which reproduces all the possi- 
ble asymptotic behaviors at early times of viable scalar-tensor cosmological 
models with a monotonic cj(<i>). This form gives an exact representation for 
most of the particular theories proposed in the literature, but also a first-order 
approximation to many other theories. In most of scalar-tensor theories, the 
comparison of our results with current observations implies very strong bounds 
on the allowed deviation from General Relativity (GR). These bounds lead 
to cosmological models which do not significantly differ from the standard 
Friedman-Robertson- Walker ones. We have found however a particular class 
of scalar-tensor theories where the expansion rate of the universe during nu- 
cleosynthesis can be very different from that found in GR, while the present 
value of the coupling function uj is high enough to ensure compatibility with 
solar-system experiments. In the framework of this class of theories, right 
primordial yields of light elements can be obtained for a baryon density range 
much wider (2.8 ^ r]io ^ 58.7) than in GR. Consequently, the usual constraint 
on the baryon contribution to the density parameter of the universe can be 
drastically relaxed (0.01 ^ Qbo ~ 1.38) by considering these gravity theories. 
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This is the first time that a scalar-tensor theory is found to be compatible both 

with primordial nucleosynthesis and solar-system experiments while implying 

cosmological models significantly different from the FRW ones. 
04.50.-Fh, 98.80.Ft, 98.80.Hw 
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I. INTRODUCTION 



From a theoretical point of view, the most natural alternatives to General Relativity (GR) 
are scalar-tensor theories which contain, in addition to the metric tensor, (/^j,, a dynamical 
scalar field, (/>, the relative importance of which is determined by an arbitrary coupling 
function c<;((/)). In recent years, this class of metric theories has received a renewed interest 
[1,2], because it provides a natural (non-fine-tuned) way to restore the original ideas of 
infiation while avoiding the cosmological difficulties coming from the vacuum-dominated 
exponential expansion obtained in GR. Scalar-tensor theories also arise in current theoretical 
attempts at deepening the connection between gravitation and the other interactions. For 
example, in modern revivals of the Kaluza-Klein theory and in supersymmetric theories 
with extra dimensions, one or several scalar fields arise in the compactification of these 
extra dimensions [3-9]. Furthermore, scalar-tensor theories may also appear as a low-energy 
limit of superstring theories [10] 

Several scalar-tensor theories have been proposed to date: 1) Brans-Dickers theory [11], 
where u = constant ^ -3/2, 2) Dime's theory [12], with u = -3/2, 3) Barker's theory 
[13], where the gravitational coupling constant is effectively constant, 4) Bekenstein's theory 
[14], with variable rest mass, and 5) Sehmidt-Greiner-Heinz-MuUer's theory [15], which 
also includes a possible mass term for the scalar field. They are all particular cases of the 
general scalar-tensor theory by Bergmann, Wagoner and Nordtvedt [16-18]. In addition, the 
seale-eovariant theory of Canuto [19,20] has a mathematical representation similar in many 
aspects to Dirac's gravitation but with a non-dynamical scalar function fi. 

The viability of a given alternative gravity theory can be analyzed by means of two 
kinds of tests [21]: those which examine its weak field limit and those which prove its full 
exact formulation. The first mainly consists of comparing the theory predictions in the 
limit of weak gravitational fields and slow motions with Post-Newtonian experiments. The 
only metric theory which is discarded by these experiments is Dirac's gravitation [21]. On 
the other hand, strong field tests consist of matching up the exact theory predictions to 
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experiments. This is mainly achieved by means of cosmological models, whose predictions 
have to be consistent with any present observation of cosmological interest. 

The astronomical data leading to the strongest bounds on the alternative gravity theories 
are the light element abundance observations, which have to be explained as an outcome of 
the primordial nucleosynthesis process (PNP). Some authors [22-26] have got a first insight 
into these PNP constraints by including a constant speed-up factor ^ in the usual GR 
expression of the universe expansion rate. Such a simple approach implies rather stringent 
bounds on ^. For example. Barrow [23] obtained 0.8 < < 1.2 while, using more recent 
reaction rates and observational data, Casas et al. [26] found an upper limit of 1.02. 

Since, in a general scalar-tensor theory, ^ is not necessarily constant at early times, 
the above approach only explores a very limited range of models. In principle, if ^ varies 
during nucleosynthesis, the resulting bounds could be very different from those derived from 
the constant-(^ approximation. Other authors [27-35] have thus preferred a more rigorous 
approach which consists of solving, from numerical computations, both the cosmological 
and nuclear equations. Since this approach follows the time evolution of ^{(j)) and all the 
dynamical functions, it has the additional advantage of determining how the PNP bounds 
result in present limits on the parameters. Using this approach, cosmological models and 
their corresponding light-element production have been analyzed in all the particular scalar- 
tensor theories proposed in the literature. The resulting bounds were always very stringent, 
implying that the only viable models are those whose predictions do not significantly differ 
from the standard GR ones up to at least temperatures of 10^° K. 

In a previous paper [36] we found all the possible early behaviors of scalar-tensor cosmo- 
logical models with a monotonic, but arbitrary, ^{(j)) function. These behaviors exhibited a 
variety of models much wider than that contained in all the previous particular cases. The 
aim of this paper is to perform a detailed numerical study of primordial nucleosynthesis 
bounds on theories presenting each one of these possible early behaviors. Such a study will 
allow us to elucidate whether the strong bounds previously found in some particular theories 
are also expected for any other one or, on the contrary, whether there exist some cases in 
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which such bounds can be considerably relaxed. To that end, we will consider throughout 
this paper a standard scenario (i.e., a homogeneous and isotropic universe with vanish- 
ing cosmological constant, without exotic particles, etc.) in the framework of scalar-tensor 
theories. 

The paper is arranged as follows. We begin outlining in Sect. II scalar-tensor theories 
and the basic equations to analyze the light-element production. We also introduce in 
that section a representation for the coupling function. Predicted abundances for a sweep of 
initial conditions and the constraints obtained from comparison with observations are shown 
in Sect. III. Finally, conclusions and a summary of our results are given in Sect. IV. 



II. SCALAR-TENSOR THEORIES 

A. Field equations and Cosmological Models 

The most general action describing a massless scalar-tensor theory of gravitation is 
[16-18] 

^=-^ li^'^- ^</',,</''")y^^^'^ + Sm (1) 
Idtt j (p 

where TZ is the curvature scalar of the metric g^^^ g = det[g^^)^ (j) is the scalar field, and ^{(j)) 
is an arbitrary coupling function determining the relative importance of the scalar field. 
The variation of Eq. (1) with respect to g^^ and (j) leads to the field equations: 

-^{4>,,;. - □</>) (2) 

{3 + 2uj)a(f> = STT(f>-Uj'(f>^a(f>'") (3) 

where uj' denotes dio/dcj) and □(/> = g^^ (f^^fi-v In addition to Eqs. (2) and (3), we have the 
standard conservation law T.^^ = 0, where T^^ is the energy-momentum tensor. 

In order to build up scalar-tensor cosmological models, we consider a homogeneous and 
isotropic universe. The line-element has then the Robertson- Walker form: 



ds' = -df + R^{t)[-^^^ + r'd^'] (4) 
i — 71 

and the energy-momentum tensor corresponds to that of a perfect fluid 

T^'- = {pj,Pl^)u^u, + Pg,, (5) 

where K = 0, ±1, R{t) is the scale factor, p and P are the energy-mass density and pressure, 
respectively, and is the 4-velocity of the fluid. The field equations (2)-(3) can be then 
written, in terms of H = R/ R and D = </>/</), as 



(3 + 2uj) 



D' (6b) 



where dots mean time derivatives. 

By eliminating R from the time-time and space-space components of Eq. (2), we obtain 
moreover the algebraic expression 

c^K Stto 9 ^ „9 „ 

R^ 3(1) 6 ^ ' 

Finally, by assuming a standard particle content, the state equation is given by the usual 
form [37] and the energy-momentum conservation law gives the usual expression for the time 
evolution of the temperature T 

dt dpi/dT 
dT " ~3H{pi + Pi/c^) 

where pi = Pb + Pe + P-y^ Pi = -Pe + P-y and subscripts 6, e and 7 refer to baryon, electron- 
positron, and photon, respectively. 

The differential equation system (Eqs. [6], [8], and the definitions of H and D) can 
be integrated by taking as variable the temperature and an explicit form for the coupling 
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function ^{(j)) (see Sect. II C). Thus, i?, (/>, i7, D and t are considered as independent 
functions . 



B. Boundary Conditions 

In order to build up realistic cosmological models, the boundary conditions are fixed, as in 
Refs. [27-30,35], by the present values of dynamical functions. By present values (hereafter 
denoted by a sub-index 0) of dynamical functions we mean their current observational values 
or their limits from observational data. 

Some of these values, namely To and i^o, can be directly measured from observations, 
whereas (/)o, -Do, and c^K/ can be expressed in terms of other better known parameters. 
As a matter of fact, the (j)o value can be obtained from loq provided that the coupling function 
iL!{(j)) has been specified. On the other hand, by using Eqs. (6a) and (7), -Do, and c^K/Rq 
can be replaced by qo = [ — RR/ R^)o = [ — [H^ + H)/H'^)o and po: 



Ho \ , , . .n- f FRW 3+^0 



and 

c^K Stt 9 cuo „9 „ , , 

Ro 300 D 

where J^o = ^-^(1^1^, = 4:TtGoPo/3H^ , and po = pm + pbo , with pRo = ye«{To)aT^, 

Pbo = 6.639 X 10~"^^(ro/2.7)"^?7io, (/efr(7o) is the effective number of relativistic degrees of 
freedom, and rjio is the present baryon-to-photon ratio in units of 10~^°. 

Note that, in order to have a real -Do, the present value of the deacceleration parameter 
must satisfy 

\ FRW 3 + CUo 1 , , 

which, when cuo ^ 1 and jFo ^ 1, implies qo ^ 

The numerical integration of Eqs. (6)-(8) needs to specify the set of parameters Hq^ To, 
cuo, qo and po, together with the usual constants (neutron mean-life, number of light neutrino 
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families) appearing in the state equation. A choice for the double-valued Dq solution given 
by Eq. (9) must be also specified. 

In order to avoid an excessive number of free parameters, we have taken the average 
values of the present photon temperature To = 2.73 it 0.02 (2(t) [38], and of the neutron 
mean-life, r„ = 889 it 4 s. [39]. The present Hubble parameter has been taken to be 50 
Km s~^Mpc~^ [40], and it has been assumed that the number of light neutrino families, N^^ 
is three, in accordance with the LEP and SLC results [41]. We will nevertheless discuss in 
Sect. Ill B how our results could be modified by the uncertainties in some of these quantities 
or by observational constraints different from those used here [42]. 

We have then considered as free parameters the present values of cuq, po (or ?7io), and 
go- Numerical computations show nevertheless that, for scalar-tensor (ST) theories which 
converge towards GR, the less stringent bounds are obtained when qo ~ q^^^ . A detailed 
discussion on this point for some particular cases can be found in [30,29,35]. To facilitate 
the discussion, we will just present here our results with _D^(go, •••) chosen to give the most 
conservative constraints on ST theories. 



C. The coupling function ijj{4>) 



In a previous paper [36] we showed that a convenient form for the coupling function io{(f)) 
is given by 

|3 + 2cu|= (3/A2)(x-^ + A;) (12) 

where 

1 - $ (if $ < 1) 

(13) 

$ - 1 (if $ > 1) 

and 1/2 < e < 2 in order to ensure that theories converge at present towards GR. Eq. (12) 
gives an exact representation for most of the particular scalar-tensor theories proposed in 
the literature [44] and, in addition, it contains all the possible early behaviors of any theory 
where wi^cf)) is a monotonic, but arbitrary, function of (f) 
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Our analysis in terms of the k and A parameters revealed the existence of scalar-tensor 
cosmological models with early expansion rates different from that found in the standard 
FRW scenarios [36]. This change in the expansion rate of the universe can modify the 
PNP in different ways. In the next section we will discuss the light element yields by 
distinguishing four main classes of viable scalar-tensor cosmological models (Table I). The 
first and the second class correspond to singular models with a monotonic speed-up factor 
^(T) = H I H^^^ faster or slower, respectively, than in GR. Non-singular models or models 
with a critical temperature where 3 -|- 2a; = constitute the third class. Finally, the last 
class corresponds to models with a non-monotonic i^(r) function. 

We know that each one of these classes can be subdivided according to other criteria as, 
for instance, the !^ value at the singularity, the increasing or decreasing evolution of 
Such analysis has been performed in our previous paper [36]. However, we will see that the 
classification using the four above classes is enough for the purposes of this paper. 

III. RESULTS 

A. Primordial Nucleosynthesis Bounds 

Light element production in scalar-tensor theories has been computed by using the up- 
dated reaction rates given by Caughlan and Fowler [45] and Smith et al. [42]. Boundary 
conditions were set as explained in Sect. II B. We will discuss here the yields of light ele- 
ments as a function of the baryon-to-photon ratio in units of 10~^°, ?7io, the present value 
of the coupling function, cuq, and the parameters A;, and A, characterizing the form of 
(see Fq. [12]). Although we now consider e = 1, we will show in the next subsection that 
our conclusions just depend very weakly on this parameter. Only the "^He, (D/H) and (D -|- 
■^IIe)/II primordial yields have been used to constraint each scalar-tensor theory. We have 
also computed the ''Li/H primordial abundance (briefiy discussed in section IIIB) but, due 
to its well-known uncertainties, it has not been used in our analysis. 
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1. Singular models with a monotonic ^{T) faster than in GR (Class-1) 



Some scalar-tensor models (see Figure 1) imply that, during all the PNP, the expansion 
rate of the universe is faster than in GR. The range of parameters defining these models can 
be seen in Table 1. Most of the particular scalar-tensor theories proposed in the literature 
are included within this category. 

Our computations show that all these theories have qualitatively similar implications on 
the primordial abundance of light elements. Before a general discussion, we first present as 
an example our results for the theory defined by = 3/2, A; = and e = 1. 

Figures 2 show the primordial abundance of "^He (denoted as 1^), (D/H), and (D -|- 
■^He)/H as a function of Uq and rjiQ. We see from these figures that, for any fixed value 
of ?7io, these abundances increase as Uq decreases. This behavior can be interpreted in the 
following way. When Uq decreases, the expansion rate !^ during nucleosynthesis is faster 
(Figure 1). The temperature at which the n/p ratio freezes out is then higher ^, and this 
in turn implies a higher light element production [46]. On the other hand, the rfw dependence 
of primordial abundances for a given value of Uq is such that the "^He production increases 
with rjio while those of D/H and (D-|-"^He)/H decrease. The reason of this behavior is that 
reaction rates increase with rfw and, consequently, the "^He production from the D and "^He 
burning is more efficient for larger values of rjiQ. 

In order to analyze the compatibility between observed and predicted light-element 
yields, we have depicted the observational bounds [42] by a thick line on the abundance 
axes of Figures 2. As can be seen from these figures, the constraints on this theory are 
essentially imposed by "^He, and (D -|- ■^He)/H. To have a right abundance of these elements, 
the conditions Uq ^ 10^^ and rfw ~ 3 are required. In particular, for very large Uq values. 



^The freezing out temperature, T*, is roughly determined by the condition: \{T^)t{T^) ~ 1, where 
A(r*) is the rate for the weak interactions changing protons into neutrons and vice versa, and t{T^) 
is the age of the universe. 
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"^He requires 1 ^ rjio ~ 3.7, while (D + ■^He)/H imposes rjio ~ 2.8. Consequently, the interval 
2.8 ^ rjio ^ 3.7 is required to have simultaneously right primordial yields for all the light 
elements. This rjio interval is even more narrow for smaller values of tOo until, for tOo ^ 10^^ 
(or, equivalently, ^iq = (^(10^°K) ^ 1.024), simultaneously right abundances are not further 
obtained for any value of rjiQ. 

Primordial nucleosynthesis bounds on other class- 1 scalar-tensor theories, with different 
A, k and e parameters, can be obtained by performing a similar computation as before. 
However, since the predicted abundances in this class of cosmological models essentially 
depend on the value (or, equivalently, on the ^iq value), the upper bound on ^iq will 
remain nearly the same whatever the theory is 

60 ^ 1.024. (14) 

On the contrary, since such a value is obtained for loq values depending on the theory 
under consideration, the loq bounds will be in general different. The k and A dependence of 
these (jJq bounds can be deduced from Figs. 3a and 3b, where we show the results for that 
value of rjio (=3) leading to the less restrictive constraints. We see from these figures that 
the larger A and the smaller A;, the higher bounds on loq are found. Such a dependence does 
not however imply small loq bounds in the limits of small A or large k values because, in 
these limits, the constraint on loq does not further decrease and it remains equal to 

tuo>102i. (15) 

Clearly, these bounds on loq are much more stringent than those obtained from the Post- 
Newtonian experiments [tOo ^ 500). They imply that any class-1 scalar-tensor theory is 
indistinguishable from GR from the beginning of PN up to the present. The allowed range 
for rjio is then similar to the usual one (2.8 ^ rjio ^ 3.7), what in turn implies the usual 
constraints on the baryonic contribution to the density parameter 

0.01 < Ofco ^ 0.09 (16) 

where the relation 
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n^oh' = 3.77 X 10-3?7io(ro/2.73)V/($o) (17) 

has been used, and where we have adopted the more conservative uncertainty 0.4 ^ h ^ 1 
[42] for the Hubble parameter (in units of 100 km s~^ Mpc~^). The correction factor /($o) 
appearing in Eq. (17) takes into account that, in scalar-tensor theories, the density param- 
eter has an additional dependence on $. However, since /($o) = $o[l — (^o/6)(-Do/-ffo)^ + 
(Do/Ho)], where (Do/Ho) — I/cuq and | $o |— ^/'^o^^ (see Eqs. 7, 9 and 12), the strong 
PNP bounds on loq (15) imply that /($o) is extremely close to unity and, consequently, its 
influence on O^o is absolutely negligible. 

2. Singular models with a monotonic ^{T) slower than in GR (Class-2) 

In models where !^ is, during nucleosynthesis, effectively slower than in GR (see Figure 
1 and Table 1 for the range of parameters A and k implying this kind of behaviors), the 
freezing-out temperature is now smaller than in GR. The smaller cuq, the smaller the 
value is, and also the light-element production. Such a Uq dependence is shown in Fig. 4 for 
the theory defined by = 3/4, A; = (cuq < and $o > !)• On the other hand, concerning 
the rjio dependence of primordial abundances, we have again that the "^He production from 
the D and "^He burning increases for larger and larger values of rjiQ. If the expansion rate is 
slower than in GR, this burning works over a longer period. Consequently, the smaller (^lo (or 
equivalently Uq for this class of models), the stronger the D/H and (D-|-"^He)/H decreasing 
with rjio is. Compatibility with observations requires again rather stringent bounds on cuq, 
(^10 and rjiQ. When Uq ^ 10^^ (or (^lo — 1), the lower and upper bounds on rfw are respectively 
imposed by (D-|-"^He)/H and "^He, like in class-1 theories. However, when cuq ~ 10^^ ((^lo < 1), 
the resulting underproduction of "^He can be balanced by increasing rjiQ. In this last case, 
the upper bound on rfw is instead imposed by D/H and, hence, this bound is expected to 
have a value somewhat larger than in class- 1 models. 

Primordial nucleosynthesis bounds on other class-2 scalar-tensor theories have been de- 
duced from similar computations. Figure 3c illustrates the case rfw = 3. When A^ is very 
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close to 1 , we get 

tuo>102i, (18) 

and the expansion rate of the universe can deviate from that obtained in GR by at most 4% 

0.96 < eio < 1. (19) 

These bounds are stronger for theories defined by smaller A values (see Fig. 3c), contrary 
to what we found in models of class-1. A k dependence is now meaningless because, in 
order to have the kind of behavior defining this class of theories, a A; = value is required. 
Finally, the allowed range of rjio is somewhat wider than that obtained in FRW cosmologies 
(2.8 ^ rjio ^ 8.4). This in turn also implies a slightly wider range for O^o 

0.01 < Ofco < 0.2 (20) 



3. Models with a cut-off: nonsingular or 3 + 2cj = at a finite temperature (Class-3) 

Nonsingular models have a maximum value of the universe temperature, T^ax^ which 
corresponds to the minimum of the scale factor. Similarly, models where 3 + 2a; vanishes at 
some finite temperature have also a maximum value of T. The smaller cuq, the smaller the 
Tmax value is. Consequently, in this class of models, the bounds on ^iq and loq are mainly 
imposed by the condition that T^ax must be high enough to allow the PNP itself. Figures 
5 show, as a function of k and A, the minimum loq satisfying this last condition. As can 
be seen from these figures, the occurrence of PN is only possible for very large values of loq 
implying an expansion rate during nucleosynthesis extremely close to that obtained in FRW 
models. Thus, in spite of its very different behavior at early times, this class of theories 
must be extremely close to GR from the beginning of PN up to the present. The required 
bounds on rjio and O^o are then very close to the standard ones (see Sect. Ill A 1). 

The A and k dependence of iOo^min is different for different intervals of these parameters. 
Theories defined by cuq > 0, $o < 1 (Fig. 5a) imply larger iOo^min values for larger A or smaller 
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k values (note that we do not write | A; |). The same dependence on A but the opposite on 
k is however found in theories defined by cuq < 0, $o < 1 (Fig. 5b). Finally, theories with 
cuo < 0, $0 > 1 (Fig. 5c) imply larger values of iOo^min for smaller k values, but they have a 
non-monotonic dependence on A. No class-3 models exist for cuq > 0, $o > 1. 

It is important to note that an appropriate choice of A (for instance, the limit A^ — > 
in Figs. 5a or 5b) and k could allow for not very large values of iOo^min (< 10^^). However, 
in that case, the bounds on loq and ^iq are those found for models with a monotonic ^{T) 
(Sections IIIAl and III A 2). 

4- Singular models with a non-monotonic ^{T) function (Class-4) 

Theories with a non-monotonic evolution of the speed-up factor ^ (see Fig. 1 and Table 
1) have an initial phase where the expansion of the universe is slower than in GR but, 
afterwards, it becomes faster than in the standard cosmology. Figs. 7 show the light 
element yields as a function of loq and rjio in the theory defined by A^ = 3/2 and k = 1/2. 
For extremely large loq values, the maximum of ^ is reached much before the PNP (Fig. 
1). Consequently the speed-up factor is greater than unity during nucleosynthesis and the 
resulting bounds on cuq, ^lo and rjio are then equal to those found for theories of class-1 (Sect. 
Ill A 1). However, we can also see from these figures that each element has another allowed 
interval for not so large values of loq. In this last interval, simultaneously right primordial 
abundances can be obtained for large rjio values. For instance, if cuq = 6 X 10^'', "^He requires 
10.0 ^ rjio ~ 15.2; D/H requires rjio ~ 11.6 and (D-|-"^He)/H imposes rjio ~ 4, what leaves 
the allowed range 10.0 ^ rjio ^ 11.6. Other values of loq in this particular theory, imply 
different bounds on rjio making possible any value between 2.8 and 12.4. The last upper 
bound on rjio needs a present value of the coupling function tOo ^ 10^^ and, hence, the theory 
can differ from GR much more than in the previous classes {^max ~ 1.6). 

The above example shows that the right primordial abundances can be obtained in the 
framework of a class-4 theory clearly distinguishable from GR. Furthermore, this class of 
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theories allows for quite large rjio values. The achievement of such rjio values can be explained 
in the following way. In the previous particular theory, the expansion rate of the universe 
at the beginning of nucleosynthesis is slower than in GR ((^lo = 0.72 if cuq = 6 X 10^'' 
and rjio = 11). Similar to what we found for Class-2 theories, there exists a tendency 
to the underproduction of "^He, which can be balanced by considering larger rjio values. 
Consequently, the upper bound on rjio is here constrained by the D/H abundance. However, 
unlike class-2 theories, the expansion rate of the universe becomes during the PNP faster 
than in GR. The D burning is not then very effective because it occurs in a shorter time 
and, hence, large rjio values are allowed. 

In order to obtain the largest rjio value allowed in the framework of this class of models, 
we must then analyze those theories with the largest ^ value during the PNP, but with 
(^10 < 1 at the beginning of PN. These two conditions are better obtained if the slope of 
(^(r), in its increasing interval, is as high as possible. Figures 8 and 9 show the A and 
k dependence of primordial abundances in this class of theories with rjio = 3. We see 
that the form of these curves is similar to that of ^{T)^ (especially Yp because it has the 
strongest dependence on ^lo). Consequently, any discussion about the behavior of ^{T) will 
be performed in terms of, for example, Yp[ijJo). The sharpness of the maximum of ^ depends 
mainly on k. The smaller A;, the sharper the maximum of ^ is. On the other hand, the 
^max value increases as A decreases while the k dependence of ^max is non-monotonic, with 
higher ^max values in the limits of very large and very small k values. In order to have a 
very large ^ value during nucleosynthesis we then need to analyze the limit of very small 
A and k values. We have performed such an analysis and we have found that, in fact, the 
allowed O^o value is considerably increased with respect to that obtained in other classes of 
scalar-tensor theories. However, the lower bound on loq is smaller as the allowed O^o is larger. 
The upper bound on the baryon contribution to the density parameter is then determined 
by the compatibility with solar-system experiments, which requires loq > 500. This is found 
for the theory A^ = 0.08, and k = 2 X 10"'', which leads to right primordial abundances with 
Vio ^ 58.7, that is 
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0.01 < ^ 1.38 



(21) 



The above result does not signify that O^o values equal or greater than unity are only possible 
for rather low values of loq. For instance, for = 0.1 and k = 10~^, (see Fig. 10), we can 
in fact always find some interval of tOo greater than 10'' allowing nevertheless for rjio values 
in the range 2.8 ^ rjio ^ 56.5. The resulting permitted interval for O^o is 



Closure of the universe by baryons is then possible in the framework of this class of scalar- 
tensor theories even when a standard composition of the universe is used [47]. The loq 
values needed to allow for O^o ^ 1 are very large and, consequently, these of theories are at 
present close enough to GR to ensure compatibility with post-Newtonian experiments. On 
the contrary, they are very different from GR at early times and imply expansion rates of 
the universe which can be several hundred of times faster (when ^ = i^maa;)than in the FRW 
cosmology. 



In the above analysis we have taken = 3, r„ = 889 s. To = 2.73, Hq = 50, and e = 1 
to avoid an excessive number of free parameters. We will now analyze how the uncertainties 
on these parameters could modify our bounds on ST theories. 

The width of the Z° particle is directly related to the number of light particles which 
can couple to the The measurements of the Z° width from results of the LFP and SLC 
e"'"e~ colliders [41] imply that the number of light neutrino species is < 3.25 at the 2cr 
level. Thus, the choice of equal to the three known species of light neutrinos does not 
seem at present a source of uncertainty in PN results [48]. 

Determinations of the neutron half-life based on the storage of ultra-cold neutrons [39] 
have considerably reduced the uncertainty in this parameter. By combining these with 



0.01 < Ofco ^ 1.33 




B. Uncertainties on the Primordial Nucleosynthesis bounds 



1. Influence of parameters 
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other recent results, Smith et al. [42] infer r„ = 889.8 it 3.6 s. We have performed some 
computations where the neutron mean-life was varied within the above allowed range. We 
obtained that, when r„ is changed from 886.2 to 893.4 s, primordial abundances vary as 
AYp = +0.0015, and AD/H = A(D+3He)/H = +0.018 xlO'^ This implies an uncertainty 
of ±0.01 for our bounds on ^iq and possible differences in a factor of two in our constraints 
on I |. The uncertainty in the neutron mean-life has then a very small effect on our 
conclusions of section III A implying that theories with a monotonic evolution of ^ are not 
significantly different from GR. The constraints on tOo remain very strong (|c<;o |^ 10^°) and 
the deviation from GR at the beginning of PN is at most 3.4% in class-1 theories, or 5% in 
class-2 theories. In the same way, our results concerning models with a non-monotonic ^{T) 
do not change significantly. This class of theories can deviate strongly from GR and allows 
for very large rjio values. The uncertainty in r„ can only modify the bounds on rjio given in 
section III A 4 by ±0.3. 

The uncertainty in the other usual PN parameters (To and Hq) has no significant effect 
on primordial abundances for a given ^iq value. The Hubble parameter can however produce 
an uncertainty in the loq value leading to that speed-up factor, (^lo, at the beginning of PN. 
For instance, computations with h = 1 lead to bounds on | cuq | which are smaller by a factor 
of 10 than those obtained with h = 0.5. The constraints on ^iq do not change and neither 
those on rjiQ. As before, the high loq bounds we have found imply that the uncertainty 
introduced on this constraint is not physically important. 

With respect to the e influence, our computations with 1/2 < e < 2 show that the 
resulting bounds on loq are slightly smaller if e is greater than unity (they decrease at most 
by a factor of ~ 10 when e is close to 2), but they increase very strongly for e values smaller 
than unity (for example, when e = 0.95 in a class-1 theory, the bound on loq typically increase 
by a factor of 10^, then implying tOo ^ 10^''). On the contrary, the bounds on (^lo and rjio do 
not depend signiflcantly on e. We then flnd that our constraints on scalar-tensor theories 
cannot be relaxed by the influence of this parameter more than by considering, for example, 
the effect of the uncertainty in Hq. 
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2. Consequences due to other observational constraints on primordial abundances 



Concerning the observational constraints on primordial abundances, we have used in 
this paper the most widely accepted ones. Recent measures of primordial abundance of 
deuterium [49] from the absorption spectra of high-redshift quasars seem however indicate a 
rather high observational lower bound for D/H (~ 1.9 X 10"'*). These measures are at present 
preliminary and cannot be reconciled with the upper bound of (D+"^He)/H (~ 9 X 10~^). 
If we however accept such a constraint for D/H and we then ignore that for (D+"^He)/H, 
the upper bound rfw is considerably reduced in GR and also in scalar-tensor theories. In 
particular, in GR and scalar-tensor theories with a monotonic i{T), compatibility with 
such observations would require rfw ^1.6 and O^o ~ 0.04, which leaves few room for non- 
luminous baryonic matter in the Universe. This effect would be however much less dramatic 
for scalar-tensor theories with a non-monotonic i{T). For example, in the theory defined 
by = 0.1 and k = 10~^ (see Eq. 22), simultaneously right abundances of "^He and D/H 
can be obtained for rjiQ ^ 12.5 and cuq ^ 5 X 10'', which implies O^o ~ 0.3. This last bound 
can be increased to O^o ~ 0.4 by considering other class-4 theories. These values of O^o are 
compatible with most of measures of the total density parameter O^o- Consequently, even 
when the high-redshift quasar measures are considered, class-4 theories allow for large rfw 
values and a universe content largely dominated by baryons. 

Finally, we remark that primordial productions of ''Li/H were computed in all runs 
of our nucleosynthesis code. However, due to its well known uncertainties, we have not 
used that element to constraint scalar-tensor theories. We note however that all the above 
constraints imply a ''Li/H production compatible with the more conservative observational 
constraint for this element (''Li/H< 1.3 X 10~^). Consequently, our bounds would not be 
modified by considering this abundance. This is also true for class-1 to class-3 theories when 
a more severe observational range is used (1.0 X 10~*° <''Li/H< 2.3 X 10~*°). However, the 
Ofco upper-bound obtained in class-4 theories would be reduced by a factor of ~ 0.5, then 
implying O^o ~ 0.65. 



18 



IV. CONCLUSIONS AND GENERAL DISCUSSION 



Primordial production of light elements has been calculated in the framework of Scalar- 
Tensor cosmological models for a sweep of initial conditions compatible with astronomical 
data. Our results and conclusions are absolutely different for theories with a monotonic 
evolution of the speed-up factor ^{T) and those with a non-monotonic ^{T). 

In the first case, our results imply that the expansion rate of the universe at the beginning 
of primordial nucleosynthesis can differ from that obtained in the usual FRW model by at 
most 4% 

0.96 < 60 ^ 1.024 (23) 

Consequently, the present value of the coupling function must be greater than 

cuo > 10'° (24) 

and, if > 1, the allowed range for the baryon density parameter is essentially the same as 
in GR 



(class-1) (25) 



(class-2) (26) 



2.8 <?7io<3.7 
0.01 < Ofco ^ 0.09 
or, if < 1, just slightly wider 

2.8 <?7io<8.4 
0.01 < 060^0.2 

These bounds, together with the cosmological evolution of Scalar- Tensor theories studied 
in paper I [36], imply that all these cosmological models are indistinguishable from the 
standard FRW ones from the beginning of PN up to the present (except for slight differences 
on the rjio bounds in class-2 models). The uncertainties on the usual input PN parameters 
can just slightly relax the above bounds and, hence, they do not modify our conclusions. 

Primordial nucleosynthesis is then a very strong test for the viability of this kind of 
gravitational theories. Very small deviations from GR during the early stages of the evolution 
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of the universe imply a light element production inconsistent with present observations. 
However, we cannot assure that these theories are also indistinguishable from GR at very 
early epochs before nucleosynthesis or with a non-standard composition of the universe 
[36,50]. 

We would like also to note that the general PN constraints given by Eq. (24) for Scalar- 
Tensor theories with a monotonic ^{T) are of the same order of magnitude as those obtained 
by us and by other authors in the framework of the particular cases proposed in the literature 
[29-35]. These PN constraints are instead much higher than those implied by the post- 
Newtonian experiments [21], and also than those found by Damour and Nordtvedt [51] 
(who estimated the matter-dominated evolution of ^{(j)) in flat scalar-tensor theories by 
assuming ^{(j)) ~ 1 at the beginning of that era). Primordial Nucleosynthesis then arises 
as a test which imposes constraints on the early value of ^{(j)) which are roughly similar 
to those found from post-Newtonian experiments for the present loq value. However, since 
iL!{(j)) increases during the universe evolution, such constraints are much stronger at the 
beginning of the matter-dominated era, and even much stronger at present. Note also that, 
by construction, other cosmological tests (measured values of i^o, 9o, To, etc. [40]) are also 
satisfied in all the models analyzed in this paper. 

Concerning theories with a non-monotonic ^{T)^ our analysis leads to very different 
conclusions. This class of theories can reproduce the right yields of light elements even 
when they are very different from GR during the radiation-dominated era. In some theories, 
the constraint on loq is similar to that imposed by solar system experiments 

uo > 500. (27) 

In such cases, the maximum value of ^ is much higher than unity 

^max ~ 350 (28) 

and the allowed range for the baryon density parameter is much wider than in GR 

2.8 < ?7io < 58.7 
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(29) 

0.01 < Ofco ^ 1.38 

Allowed Ofco values greater than unity can also be obtained in other class-4 theories 
implying a much larger present value of the coupling function loq. For instance, theories as 
that of Eq. (22), imply tOo ^ 10'' while the upper limit of O^o is only slightly smaller than 
that given by Eq. (29): O^o ~ 1.33. The universe closure by baryons is then possible in the 
framework of these theories without requiring a non-standard composition for the cosmic 
gas [47]. 

Eurthermore, we have shown in Sect IIIB that, if we accept some recent measures of 
D/H from high-redshift quasars, the resulting upper bound on O^o would be drastically 
reduced in the framework of GR as well as in theories with a monotonic (,{T): O^o ~ 0.04. 
On the contrary, the O^o bounds in models with a non-monotonic ^{T) would be much less 
severely affected by such measures. They would be still much larger (O^o ~ 0.4) than the 
usual one. On the other hand, it is also important to note that all the bounds given by 
the above expressions also imply a right abundance of ''Li/H when the most conservative 
observational constraints (''Li/H< 1.3 X 10~^) are considered. Nevertheless, if we use more 
severe observational constraints (1.0 X 10~^° ^ ''Li/H ^ 2.3 X 10~^°), the upper bound on 
Ofco imposed by the class-4 theories would be O^o ~ 0.65, which is still considerably wider 
than the usual one. 

This is the first time that a class of scalar tensor theories is found to be compatible 
both with primordial abundances and with solar-system experiments while remaining clearly 
distinguishable from GR. This last feature is not only manifested in the much wider allowed 
range for rjio but also in the cosmological evolution of the universe during the PNP and, in 
the extreme case of Eq. (21), during all the radiation dominated era. 
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FIGURES 

FIG. 1. Expansion rate ^ = H/H^^^ as a function of T for different types of cosmological 
models. Solid lines correspond to ujq values smaller than those of dashed-lines. 

FIG. 2. Primordial abundances of a) "^He, b) D/H, and c) (D+"^He) /H, as a function of loq. The 
theory shown in this figure is defined by cjq > 0, > l,k = 0, = 3/2. Solid line corresponds to 
r]iQ = 1, dashed line to 7710 = 3, and dotted line to 7710 = 5. 

FIG. 3. A, k and ujq dependence of Yp in theories with a monotonic ^{T). Theories shown in 
this figure are defined by a) cjq > 0, > 1, A; = 3 and A^ = 1.5 (solid lines), 0.75 (dashed lines), 
and 0.01 (dotted lines); b) cjq > 0, $0 > 1, A^ = 3/2 and k = 26.5 (solid lines), 10 (dashed lines), 
and 3 (dotted lines); and c) cjq < 0, > 1, A; = and A^ = 0.75 (solid lines), 0.5 (dashed lines), 
and 0.3 (dotted lines). 

FIG. 4. Primordial abundances of a) "^He, b) D/H, and c) (D+"^He) /H, as a function of ujq. The 
theory shown in this figure is defined by cjq < 0, > l,k = 0, A^ = 3/4. Solid line corresponds to 
r]iQ = 1, dashed line to 7710 = 3, and dotted line to 7710 = 5. 

FIG. 5. A dependence of the smallest value of ujq needed to have a high enough temperature 
for the PNP in theories with a cut-off. a) ujq > 0, < 1, and k = 0.5 (solid line), -2.5 (dashed 
line), and -23.5 (dotted line), b) ujq < 0, < 1, and A; = 10 (solid line), 1 (dashed line), and -8 
(dotted line), c) ujq < 0, > 1, and k = -0.5 (solid line), -3 (dashed line), and -5 (dotted line). 

FIG. 6. Primordial abundances of a) "^He, b) D/H, and c) (D-|-"^He) /H, as a function of loq. The 
theory shown in this figure is defined by ujq > 0,^o > l,k = 1/2, A^ = 3/2. Solid line corresponds 
to 7710 = 3, dashed line to 7710 = 5, and dotted line to 7710 = 10. 

FIG. 7. A, k and ujq dependence of Yp in theories with a non-monotonic ^{T). Theories shown 
in this figure are defined by a) ujq > 0,^o > l,k = 0.01 and A^ = 1. (solid lines), 0.4 (dashed lines), 
and 0.3 (dotted lines); and b) loq > 0, > 1, A^ = 3/2 and k = 1. (solid lines), 1/2 (dashed lines), 
and 0.01 (dotted lines). 
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FIG. 8. Primordial abundances of a) "^He, b) D/H, and c) (D+"^He)/H, as a function of 7710. 
The theory shown in this figure is defined by loq > 0, $0 > 1, = 10~^ and k = 10~^. Lines 
correspond to loq = 10'' (solid lines), 5 X 10'' (dashed lines), and 10^^ (dotted lines). The last case 
gives cosmological models close to the FRW ones 
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